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ABSTRACT 



Objective: Autism Spectrum Disorder (ASD) has been associated with a complex pattern of increases and 
decreases in resting-state functional connectivity. The developmental disconnection hypothesis of ASD poses 
that shorter connections become overly well established with development in this disorder, at the cost of 
long-range connections. Here, we investigated resting-state connectivity in relatively young boys with ASD 
and typically developing children. We hypothesized that ASD would be associated with reduced connectivity 
between networks, and increased connectivity within networks, reflecting poorer integration and segregation 
of functional networks in ASD. 

Methods: We acquired resting-state fMRI from 27 boys with ASD and 29 age- and IQ-matched typically 
developing boys between 6 and 16 years of age. Functional connectivity networks were identified using 
Independent Component Analysis (ICA). Group comparisons were conducted using permutation testing, 
with and without voxel-wise correction for grey matter density. 

Results: We found no between-group differences in within-network connectivity. However, we did find 
reduced functional connectivity between two higher-order cognitive networks in ASD. Furthermore, we 
found an interaction effect with age in the DMN: insula connectivity increased with age in ASD, whereas it 
decreased in typically developing children. 

Conclusions: These results show subtle changes in between network connectivity in relatively young boys 
with ASD. However, the global architecture of resting-state networks appeared to be intact. This argues against 
recent suggestions that changes in connectivity in ASD may be the most prominent during development. 

© 2014 The Authors. Published by Elsevier Inc. 
This is an open access article under the CC BY-NC-ND license 
(http://creativecommons.Org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Autism Spectrum Disorder (ASD) is a developmental disorder char- 
acterized by restricted and repetitive behaviour, as well as impaired 
social interaction and communication (American Psychiatric Associa- 
tion, 1994). ASD has been associated with changes in brain connectiv- 
ity (Belmonte et al, 2004; Courchesne and Pierce, 2005; Geschwind 
and Levitt, 2007; Minshew and Keller, 2010), where some connec- 
tions are reported to be weaker (Just et al, 2004; Just et al, 2007; 
Kana et al, 2007; Kleinhans et al, 2008; Villalobos et al, 2012; von 
dem Hagen et al, 2013) and others stronger (Koshino et al, 2005; 
Mizuno et al, 2006; Turner et al, 2006; Shih et al, 2010; Di Martino 
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et al, 201 1 ) in individuals with ASD compared to typically develop- 
ing controls. The developmental disconnection model of ASD suggests 
that changes in functional connectivity in individuals with ASD fol- 
low a pattern of local over-connectivity and long-range (Belmonte et 
al, 2004; Courchesne and Pierce, 2005; Geschwind and Levitt, 2007), 
or even global under-connectivity (Just et al, 2004; Just et al, 2007). 
This aberrant pattern of connectivity may then reflect deficits in the 
fine-tuning of networks during development in ASD, a process that 
is typically characterized by increased integration within networks 
and segregation between networks (Fair et al, 2007; Fair et al, 2009; 
Supekar et al, 2009). It has been suggested that development of func- 
tional connectivity in ASD may be characterized by atypical rather 
than by delayed integration or segregation of functional networks 
during childhood (Di Martino et al, 2011; Washington et al, 2014), 
where the developmental trajectories may differ even between func- 
tional networks (Padmanabhan et al. 2013). 

is an open access article under the CC BY-NC-ND license (http://creativecommons.org/ 
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Many studies have investigated functional connectivity during 
cognitive tasl<s (Just et al., 2004; Koshino et al, 2005; Mizuno et al, 
2006; Turner et al, 2006; Just et al, 2007; Kana et al, 2007; Kleinhans 
et al, 2008; Shih et al, 2010; Villalobos et al, 2012), although more re- 
cently investigators have also started to study functional connectivity 
during rest in ASD (Cherkassky et al, 2006; Monk et al, 2009; Assaf 
et al, 2010; Weng et al, 2010; Di Martino et al, 2011; Ebisch et al, 
201 1 ; Wiggins et al, 201 1 ; von dem Hagen et al, 2013 ; Washington et 
al, 2014). In the general population, resting state networks identified 
thus far appear to resemble cortical functional 'task' networks, such as 
the sensorimotor, visual, and executive networks (Damoiseaux et al. 
2006; Seeley et al. 2007; Vincent et al. 2008; Laird et al. 201 1 ; Shirer et 
al. 2012). The only network that does not resemble a 'task' network is 
the Default Mode Network (DMN), which includes the posterior cin- 
gulate gyrus (PCC), retrosplenial cortex, lateral parietal cortex, medial 
prefrontal cortex (mPFC), superior frontal gyrus (SFG) and temporal 
lobe (Greicius et al, 2003; Fox et al, 2005). The cerebellum seems 
to be involved in many of these intrinsic networks, including higher 
order cognitive networks and the DMN (Allen et al, 2005; Habas et 
al., 2009). 

Even though ASD is a developmental disorder, only a few stud- 
ies have investigated resting-state functional connectivity in chil- 
dren with ASD. Most have reported reduced functional connectiv- 
ity between core areas of the DMN, especially between the pos- 
terior PCC and mPFC (Di Martino et al. 2013), and as such seem 
to suggest decreases in the development of long-range functional 
networks in children with ASD (Assaf et al, 2010; Wiggins et al, 
2011;Padmanabhan et al. 2013; Washington et al, 2014). In con- 
trast, connectivity within a prefrontal DMN node was reported to 
be increased in ASD (Washington et al, 2014). Increased functional 
connectivity between different striatal areas during rest was also re- 
ported in young children with ASD (Di Martino et al., 2011), as well 
as reduced connectivity between the insula, and the somatosensory 
cortices and amygdala (Ebisch et al, 2011). These changes in mostly 
long-range functional connectivity have been associated with social 
deficits in ASD, where more severe problems in social interaction and 
communication deficits were related to lower connectivity within 
DMN areas (Assaf et al., 2010; Lynch et al., 2013). 

In the current study we investigated resting-state networks in 
a group of relatively young, high-functioning boys with ASD and a 
group of age-matched, typically developing control boys. In accor- 
dance with the developmental disconnection hypothesis, we hypoth- 
esized that there would be increased connectivity within networks, 
and decreased connectivity between networks in boys with ASD, re- 
flecting atypical development of functional networks. Furthermore, 
we hypothesized that we would find these between-group differences 
specifically in the DMN and networks related to higher-order cogni- 
tive functioning and that they would also involve the cerebellum, 
given the well-established cognitive impairments in ASD. 

2. Materials and methods 

2.1. Participants 

Forty high-functioning right-handed male children and young 
adolescents with a diagnosis within the DSM-IV autism spectrum 
were recruited through the Department of Psychiatry at the Univer- 
sity Medical Center in Utrecht and through advertising. The clinical 
diagnosis was confirmed with a research diagnosis by a qualified re- 
searcher using the Autism Diagnostic Interview — Revised (ADI-R 
(Lord et al., 1994)). The children with ASD were either medication 
naive or using psychostimulant medication (see: Table 1). No other 
forms of psychoactive medication were permitted in this study. The 
group of children with ASD that was on stimulant medication was 
instructed not to take their medication 24 h prior to the MRI-scan. 

Forty typically developing subjects were selected from the pool 



of volunteers participating in studies by our lab. All subjects were 
screened by phone interview, to confirm the absence of major neuro- 
logical or psychiatric disorders, as well as the absence of psychiatric 
conditions in their first-degree relatives. None of the control subjects 
were using any form of psychoactive medication. Parents of the typi- 
cally developing subjects participated in a semi-structured interview 
session (Diagnostic Interview Schedule for Children — Parent Version 
(DISC-P) (Costello et al, 1985)) with a trained researcher to confirm 
the absence of any psychiatric condition in the subject. 

Subjects were matched for age. Tanner stage, IQ and socio- 
economic status (see: Table 1). Tanner stage was assessed in order 
to preclude group differences in pubertal development. Full scale IQ 
was assessed for all participants using the Wechsler Intelligence Scale 
for Children (WISC-III (Kort et al, 2005)). 

2.2. fMRI acquisition 

All participants under 13 years of age participated in a mock MRI- 
scanner practice session in our lab. The purposes of this practice ses- 
sion were to acclimate the children to the scanner environment, to 
reduce potential anxiety and to increase compliance to reduce the 
risk of motion artefacts (Durston et al, 2009; Langen et al, 2009). 
Children only participated in the actual MRI-scan after a successful 
practice session. Older subjects were also offered such a session, in 
which all of them participated. 

Data were acquired using a 3.0 T Phillips Achieva MR scanner 
(Philips Medical Systems, Best, The Netherlands). A high-resolution 
Tl -weighted image was acquired to allow for spatial normalization 
and visualization (TR/TE = 10/4.6, flip angle = 8°, matrix 304 x 299, 
voxel size 0.75 x 0.75 x 0.8). During the acquisition of the Tl- 
weighted image, all children were watching a movie. Resting-state 
functional images were then collected in a single block of 294 dy- 
namics with a 2D-EPI SENSE sequence (TR/TE = 2000/35, flip an- 
gle = 70°, matrix 68 x 66, voxel size 3 x 3 x 3.5), with a total 
duration of 10 min. Participants were instructed to focus on a fixation 
cross during the procedure. 

2.3. Image preprocessing 

Resting-state fMRI data were preprocessed using SPM8 (Well- 
come Dept. of Cognitive Neurology, http://www.fil.ion.ucl.ac.uk). First, 
rigid-body realignment was performed, followed by unwarping to 
remove residual variance related to movement. Only movement not 
exceeding the total linear displacement of 3 mm or the size of 1 voxel 
over the course of the session was allowed. Hence, 13 children with 
ASD and 11 typically developing children were excluded. Average 
head motion in the group carried forward for analysis was 1.24 mm. 
Framewise displacement (as computed by Power et al, 2012) and 
Root Mean Square mean motion (as calculated by Van Dijk et al, 
2012) did not differ between groups, and did not correlate with age 
(see: Table 2). The next step was co-registration of the functional and 
anatomical images. Anatomical images were segmented into grey and 
white matter. Both the functional and anatomical images were then 
normalized to standard stereotactic space (Montreal Neurological In- 
stitute (MNI) — template). Finally images were smoothed with an 
8 mm^ FWHM Gaussian kernel, and a high-pass temporal filter was 
applied at 0.008 Hz to remove low frequency drift (cut-off = 128 s). 

2.4. Data analysis 

Group Independent Component Analysis (ICA) was performed 
(Calhoun et al., 2001) on the preprocessed functional images using 
the Group ICA of fMRI Toolbox (GIFT, http://icatb.sourcefourge.net, 
version 2.e). ICA is a model free method that decomposes functional 
imaging data into spatially independent, but temporally coherent 
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Table 1 

Demographic data for the sample. 





ASD(N = 27) 


TDC(N=29) 


P 


/Age lllCdll ^ DU ^IdllgCJ 




iZi.z, m z,. 1 y\j.'-tu.o J 


520 


Total IQmean ± SD (range) 


114.0 ± 14.4(88-145) 


113.9 ± 15.0(80-145) 


.993 


Parental education^ mean ± SD (years) 


13.6 ± 2.1 


14.3 ± 1.9 


.193 


Heighti^ mean ± SD 


153.1 ± 15.7 


157.7 ± 13.9 


.618 


Weighf^ mean ± SD 


43.1 ± 16.6 


45.9 ±11.3 


.564 


Tanner stage*^ mean ± SD 


1.8 ± 1.3 


2.3 ± 1.4 


.226 


ADl-R social mean ± SD 


19.9 ± 4.2 






ADl-R communication mean ± SD 


15.7 ± 4.2 






ADl-R repetitive mean ± SD 


4.9 ± 2.5 






Use of psychostimulant medication 


12/27 







ASD: Autism Spectrum Disorder; TDC: Typically Developing Controls; IQ: Intelligence Quotient; ADl-R: Autism Diagnostic Interview — Revised; 
"Information was not available for 1 TDC and 7 ASD participants. 
'^Information was not available for 4 ASD participants. 
•^Information was not available for 3 ASD participants. 
''Information was not available for 3 ASD participants. 



Table 2 

Motion characteristics. 





ASD(N = 27) 


TDC (N = 29) 


P 


Framewise displacement^ mean ± SD 


0.14 ± 0.05 


0.15 ± 0.09 


.550 


RMS motion!^ mean ± SD 


0.08 ± 0.05 


0.09 ± 0.03 


.467 



"Per Power et al. (2012). 
'^Pervan Dijket al.(2012). 



brain networks (Calhoun et al, 2001 ). To minimize the effects of resid- 
ual motion or any other artefacts, ICA was first run on each subject 
individually. After removing artefacts on the individual subject level 
all 56 subjects were included in a Group ICA analysis. The steps are 
described below. 



2.5. Single subject analysis 

The number of network components was estimated per subject us- 
ing minimum description length (MDL, Li et al, 2007). Independent 
components (ICs) were estimated using the Infomax algorithm (Bell 
and Sejnowski, 1995), which was repeated 20 times using ICASSO 
[http://research.ics.aalto.fi/ica/icasso) to reach an estimate of com- 
ponent reliability (Himberg et al, 2004). Finally, images were back- 
reconstructed (GICA3, Erhardt et al, 2011) and converted for visual- 
ization to reflect percent signal change. 

Following the single-subject ICA analysis, all spatial maps were 
inspected for the presence of obvious artefacts (e.g. scanner artefacts, 
edges). Furthermore, individual motion parameters were temporally 
regressed with the individual component maps to exclude compo- 
nents showing signal variation related to motion. Components were 
discarded if their correlation with the motion parameters exceeded 
a Z-transformed of 1.96. All nuisance components were then re- 
moved from the preprocessed resting-state data, resulting in cleaned 
images for all subjects. 



2.6. Group analysis 

The cleaned data of all 56 subjects was fed into a group ICA analysis 
using the Infomax algorithm, which was repeated 20 times in ICASSO 
to again ensure the reliability of the component estimation. Data from 
the whole group, both patients and control subjects, was decomposed 
into 25 independent components in order to represent robust large- 
scale networks (Beckmann et al. 2005; Greicius et al. 2007). Individual 
subject components were again back reconstructed using GICA3, and 
finally time-courses were scaled using maximum intensity value, and 
spatial maps using the standard deviation of the time-courses. 



2.7. Identification and selection of networks of interest 

Three components were discarded as they showed a spatial corre- 
lation of > 0.025 with either the white matter or the cerebrospinal 
fluid and could thus be considered as noise. Further identification 
of resting-state functional networks was performed through cross- 
correlation with established templates from Shirer et al. (2012). Com- 
ponents that represented primary sensory networks or noise were not 
carried forward to statistical analyses. Ten components were selected 
for analysis as they showed a moderate to good (>0.30) correlation 
with the templates and they represented Default Mode (ICs 8, 17, 19 
and 23), frontoparietal (ICs 22 and 24), striatal (ICs 4 and 18), cingu- 
late gyrus (IC 16) and cerebellar-thalamic (IC 10) components (see: 
Fig.l). 



2.8. Group differences in within network connectivity 

Group differences in within network functional connectivity were 
tested for the 10 selected networks using a general linear model with 
5000 permutations in Randomise (FMRIB Software Library randomise 
v2.9), with age as a covariate. Additionally, the interaction between 
diagnosis and age was tested within networks. Results were corrected 
for multiple comparisons, using Familywise Error (FWE) at p < 0.05 
and threshold-free cluster enhancement (Smith and Nichols, 2009). 
Initially, medication status was added as a covariate in the design. 
However, as this had no significant effect on the results, it was re- 
moved in the final analyses. 

All subjects' Tl -weighted grey matter segmentations were used 
to construct an average grey matter mask that was used to constrain 
the analysis to grey matter. Finally, for components with significant 
between-group effects, the analysis was repeated with voxelwise grey 
matter density as an additional covariate in the design. 



2.9. Group differences in between network connectivity 

The temporal correlations between the ten selected networks 
were computed as a measure of between network connectivity 
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Fig. 1. Spatial maps of all independent components (ICs) containing neuronal signal. ICs related to noise are not displayed. MNl coordinates are provided below each spatial map. 



Strength using FSLNets vO.3 (FMRIB Software Library). First, Pear- 
son's r correlation coefficients were calculated between the individ- 
ual average time-courses for each of the ten selected networks. The 
resulting correlation matrices were then transformed to Fisher's z- 
values. To test for group differences, the individual correlation matri- 
ces were fed into a General Linear Model with 5000 permutations in 
Randomise (FMRIB Software Library randomise v2.9), with age as a 
covariate. Next the interaction between diagnosis and age was tested 
within networks. Results were corrected for multiple comparisons, 
using Familywise Error (FWE) at p < 0.05. A schematic of functional 
network hierarchy is obtained by the clustering of components driven 
by the direction and strength of the correlations between any two of 
the components in this analysis. 



3. Results 

3. 1 . Within network connectivity 

There were no between-group differences within network bound- 
aries (see: Fig. 2). However, we found an interaction between group 
and age within the Default Mode Network: In the right insula, within 
network connectivity increased with age for the boys with ASD, 
whereas it decreased with age for the typically developing boys (IC 
19, Pfwe = 0.023; MNI-coordinates: x = 13,y = 33, z = 18; see Fig. 3). 
This effect persisted after correction for grey matter density. 



The analysis method uses a whole-brain approach and does there- 
fore not restrain the within network analyses to the boundaries of 
that network. Therefore, it is possible that clusters are detected out- 
side the network boundaries that show between group differences in 
connectivity with that same network (Voets et al, 2012). Compared 
to typically developing children, children with ASD showed reduced 
connectivity between a Default Mode component (IC 19) and a clus- 
ter in right superior frontal gyrus/paracingulate (pfwe = 0.031; MNI 
coordinates: x = 25, y = 56, z = 20; Fig. 4A), which persisted after 
correcting for grey matter density. Connectivity was also reduced for 
individuals with ASD between a frontoparietal network (IC 22) and 
a small cluster in the right supramarginal gyrus (pfwe = 0.045; MNI 
coordinates: x = 9, y = 29, z = 25; Fig. 4B). However, this effect did not 
survive FWE-correction for multiple comparisons after correction for 
grey matter density (Puncorr = 0.0002). Furthermore, the children with 
ASD showed reduced connectivity between a striatal network (IC 4), 
and a cluster in the left middle temporal gyrus (MTG) (pfwe = 0.018). 
This effect persisted after correction for grey matter density. 

3.2. Between network connectivity 

The strength of connectivity between a frontoparietal network (IC 
22) and a cingulate gyrus network (IC 16) was significantly lower 
for children with ASD in comparison to typically developing children 
(Pfwe = 0.042 ; see Fig. 5 ). Within the group of children with ASD these 
two networks showed no significant relation, whereas in typically 
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Fig. 2. Independent components are similar for ASD and typically developing controls (TDC). Separate mean spatial maps for both diagnostic groups (same coordinates as Fig. 1 ). 
No significant differences between groups were found in these networks. 



developing children these two networks were positively correlated. 



4. Discussion 

In this study we set out to investigate within and between net- 
work connectivity during rest in high-functioning boys with ASD and 
matched typically developing controls. Using group ICA, we identified 
ten resting-state networks related to higher-order cognitive function- 
ing. We found most major adult resting-state networks to be present 
in our group of relatively young children. We found no differences in 
within-network connectivity between boys with ASD and controls, 
although there were subtle changes in between-network connectiv- 
ity. 

Although we did not find any group differences within the net- 
works of interest, our approach permitted us to look beyond the 



boundaries of the intrinsic connectivity networks: it permitted test- 
ing all voxels in the brain that correlated with the network in order to 
detect differences in connectivity with areas that were less strongly 
connected with it. Using this method, we found subtle reductions in 
connectivity between the DMN and an area in the prefrontal cortex in 
ASD, as well as between a higher-order frontoparietal network and a 
parietal region. The first finding is consistent with earlier reports on 
reduced connectivity between superior frontal and more posterior 
regions of the brain in adults (Cherkassky et al, 2006; Monk et al., 
2009; Weng et al, 2010; Wiggins et al, 201 1 ) and children with ASD 
(Di Martino et al. 2013; Washington et al, 2014). It is noticeable that 
these areas are adjacent to the networks. It seems plausible that not 
all networks are fully mature at this young age (Fair et al, 2008; Fair 
et al, 2009; Jolles et al, 2011; de Bie et al, 2012) and while these 
areas are not as strongly connected at this young age, they may be 
incorporated into the network with development (Fair et al, 2007; 
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Fig. 3. Developmental differences in within network connectivity in insula. Panel A 
shows the DMN (IC 19: shown in heat colours), including the right insula (circled), 
where there was an interaction effect of age and diagnostic status on within network 
connectivity. Panel B shows this interaction, with beta-values from the right insula 
cluster plotted against age (orange for ASD (r^ = .353 ) and grey for typically developing 
children (r2 = .253)). 




ASD TDC 

Fig. 4. Reduced connectivity in higher order networks in ASD. Panel A shows reduced 
connectivity in ASD between the DMN (IC 19: heat colours) and a cluster in the right 
superior frontal gyrus (SFG; blue). Panel B shows the same for a frontoparietal network 
(IC22: heat colours) and a cluster in the right supramarginal gyrus (blue) (both Fami- 
ly wise Error corrected atp < 0.05). On the right, beta values for both blue clusters are 
displayed, representing the strength of connectivity for both diagnostic groups. 



Fair et al, 2008; Fair et al, 2009; Supel<ar et al., 2009; Power et al., 
2010; Supekar et al, 2010). If interpreted in this way, these results 
could be taken to suggest that in children with ASD, the process of 
network development may be altered, or delayed. Further support 
for this hypothesis comes from our finding of reduced connectivity 
between networks, specifically between a frontoparietal network and 
cingulate network. This ties in to other reports on deficits in cogni- 
tive control and lower activity in the inferior parietal lobe and dorsal 
ACC in adults with ASD during performance in cognitive control tasks 



DMN DMN BG BG DMN FP FP DMN Cb ACC 

ICI7 IC23 IC4 ICI8 ICS IC22 IC24 ICI9 ICIO ICI6 



B DMN 




Fig. 5. Limited differences in between network connectivity in ASD. Panel A shows the 
functional hierarchy of correlated and anti-correlated independent components (ICs) 
for the complete group of subjects. Each component that clusters under the blue lines 
(left) has a negative correlation with the components clustering under the red lines 
(right). Within these two groups of components (red and blue), the two components 
showing the strongest positive correlation are depicted on the bottom left. The further 
a component is placed to the right, the weaker its correlation with the other com- 
ponents. In panel B, correlations between component time-series are displayed with 
negative correlations in blue, and positive correlations in red. Correlations for the ASD 
group are shown under the diagonal, whereas correlations for the TDC group are above 
the diagonal. There is a between group difference in between-network connectivity 
between a frontoparietal network (IC 22) and the ACC network (IC 16) (demarcated 
with asterisk; p < 0.05, FWE corrected for multiple comparisons). DMN: Default Mode 
Network, BG: basal ganglia, FP: frontoparietal, Cb: cerebellum, ACC: anterior cingulate 
cortex. 



(Kana et al., 2007; Solomon et al, 2009; Agam et al, 2010). 

Further, we found a developmental change in connectivity in ASD: 
connectivity within the DMN increased with age in subjects with ASD 
in the right insula, whereas it decreased with age in typically devel- 
oping controls. This appears to be in contrast to earlier reports on 
insula function in adults with ASD showing reduced anterior insula 
activity during a variety of cognitive tasks and resting-state (Uddin 
and Menon, 2009; Ebisch et al, 2011; von dem Hagen et al, 2013). 
The insula is a functionally heterogeneous region (Seeley et al, 2007; 
Menon and Uddin, 201 0; Xue et al, 2010; Kelly et al, 2012), with clear 
functional differentiation between the anterior and posterior insula 
(Craig, 2002; Critchley et al, 2004). One explanation could be that 
increased activity in this region is compensatory, in reaction to less 
effective connectivity of the anterior insula. 

Finally, even though we report subtle regional differences in be- 
tween network connectivity in ASD, the global architecture of the 
resting-state networks seems intact in our relatively young sam- 
ple of children with ASD. These findings correspond to recently re- 
ported findings of largely typical resting-state functional connectivity 
in adults with ASD (Tyszka et al., 2013). However, our results chal- 
lenge the suggestion that differences in functional connectivity may 
be most pronounced during development (Anderson et al, 201 1 ), fol- 
lowed by normalization into adulthood (Tyszka et al, 2011; Tyszka 
et al, 2013). Using a similar data-driven approach, our data indicate 
that although there are some changes in communication between 
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networks, functional connectivity during rest is largely intact in chil- 
dren with ASD. 

There are some limitations that need to be taken into considera- 
tion. First, motion is a known confounder in fMRI research with chil- 
dren, especially in connectivity analyses (Power et al, 2012; Van Dijk 
et al, 2012). We addressed this by excluding subjects with excessive 
motion and stringently removing motion artefacts from the data at the 
single-subject level. Second, 12 out of our 27 subjects with ASD were 
on stimulant medication. In order to minimize the effects of medica- 
tion, subjects were instructed not to take their medication 24 h prior 
to scanning. This strategy appears to have been effective, as medica- 
tion status had no significant effect on the results. Third, we could not 
actively monitor subjects during the acquisition of the resting-state 
scan. However, none of the subjects reported falling asleep during 
the scan during debriefing. Finally, we acknowledge that our sample 
size is modest compared to collaborative resting-state fMRI datasets, 
such as the ABIDE-initiative. Nonetheless, the subjects included here 
represent an independent sample, and such datasets will be essen- 
tial, not only for replicating findings from such datasets, but also for 
generating new hypotheses to be tested with them. 

In conclusion, although we report subtle changes in between net- 
work connectivity in young children with ASD, the global architecture 
of resting-state networks appears to be intact. This argues against re- 
cent suggestions that changes in connectivity in ASD may be most 
prominent during development. 
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